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Abstract 

We present the one-loop matching condition for the unpolarized and polarized generalized quark 
distributions in the nonsinglet case. The matching condition links the quasi distributions defined 
in terms of spacelike correlators at finite nucleon momentum to the light cone distributions, and it 
is useful for extracting the latter from the former in a lattice QCD calculation. Our results show 
that at one-loop and leading power accuracy the matching for the light cone generalized quark 
distribution H (H) is nontrivial, whereas no matching is required for E {E). Therefore, E (E) can 
be smoothly approached by its quasi counterpart in the large momentum limit. We also present 
the matching for the distribution amplitude of the pion. 
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I. INTRODUCTION 


One of the important goals of quantum chromodynamics (QCD) is to understand the in¬ 
ternal structure of nucleons in terms of the fundamental degrees of freedom of QCD -quarks 
and gluons. The parton distribution functions (PDFs) play a crucial role in characterizing 
the nucleon structure. They are dehned as the forward hadronic matrix elements of light 
cone correlations, and they describe the momentum distributions of quarks and gluons in¬ 
side the nucleon. In recent years, their generalization to nonforward kinematics, known as 
generalized parton distributions (GPDs) also received considerable attention (for re¬ 

cent reviews on GPDs see e.g. |5|-l7|). In contrast to parton distributions, the GPDs encode 
more information about the internal structure of nucleons, and can be viewed as a hybrid 
of parton distributions, form factors and distribution amplitudes. Theyplayed an impor¬ 
tant role in providing a three-dimensional spatial picture of the nucleon |8| and in revealing 
the spin structure of the nucleon [I|. Experimentally, GPDs can be accessed in exclusive 
processes such as deeply virtual Gompton scattering or meson production. However, de¬ 
hned as nonlocal light cone correlations, they are rather difhcult to access by lattice QGD 
simulations. 

Recently, a direct approach to accessing parton distributions and related quantities has 
been proposed |i-[l9|- According to this approach, the light cone parton distribution can be 


studied by investigating the large momentum limit of a quasi parton distribution, which is 
a time-independent spacelike correlation and thus can be simulated on a Euclidean lattice. 
The light cone distribution is then recovered from the quasi one by a factorization formula or 
matching condition. This procedure in principle applies not only to parton distributions, but 
also to other quantities dehned on the light cone. In Ref. [l^, we presented a factorization 
formula connecting the light cone and quasi parton distributions and proved its validity up 
to one-loop order, where we showed that the quasi and light cone parton distributions have 
the same collinear singularities, and the matching factor connecting them is sensitive to UV 
physics only. 

In this paper, we consider the one-loop matching for GPDs. In particular, we focus on the 
unpolarized GPDs and and the polarized ones H{x,^,t) and E{x,^,t), 

which are dehned in terms of the following matrix elements 


Eg{x,^,t) = J (/l^(-|)7^^(-|,|)V'(|)b%+=o,^-x=o 

1 /\ 

[H{x, t)u{p")-f+u{p') + E{x, t)u{p") " uip ')], 


2p+ 

dz 


= ^ [H{x, t)u{f)j^'j^u{p') + E{x, t)u{f) ’ 

where L(—|, |) is the gauge link along the light cone and 


( 1 ) 






t = A\ e = 


pll+ _l_ p/+ 


( 2 ) 


In the limit ^,t ^ 0, H and H reduce to the usual unpolarized and polarized parton 
distributions, while the information encoded in E and E cannot be accessed since they are 
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multiplied by the momentum transfer A. Only in exclusive processes with a hnite momentum 
transfer can E and E be probed. 

We will study the unpolarized and polarized GPDs as well as their quasi counterparts 
dehned in terms of spacelike correlations, and we will compute the one-loop corrections. 
Based on the one-loop results, we then propose a factorization formula for quasi GPDs, 
and extract the matching factors relating them to the light cone GPDs. The matching for 
the GPD H {H) turns out to be similar to that for the parton distribution, whereas the 
matching for E (E) is trivial since, as we will show in this paper, the quasi and light cone 
dehnition yields the same result for E (E) at one-loop and leading power accuracy. This 
implies that the light cone GPD E [E) can be smoothly approached by the large momentum 
limit of its quasi counterpart; hence, its simulation on the lattice is relatively simple. As a 
related quantity, we also present the matching for the distribution amplitude of the pion. 

The rest of this paper is organized as follows. In Section II, we present the dehnitions 
of quasi GPDs and our conventions. In Section III, the results of our one-loop calculation 
for the unpolarized and polarized GPDs are given. The factorization formula for the quasi 
GPDs is presented in Section IV, where the one-loop matching factors are given. We also 
present the one-loop matching condition for the pion distribution amplitude. Section V 
contains our conclusions. 


II. QUASI GPDS AND CONVENTIONS 


The quasi GPDs are dehned in full analogy to the light cone ones, and they can be 
extracted from the following matrix elements dehned on a spacelike interval along the 2 ; 
direction T3 




dz 

—f 

dvr 










-u 


(p')]. 








= ^ t, P^)m(p")7^7^w(p') + £(x, t, p")u{p'') • 


( 3 ) 


The gauge link L points along the 2 ; direction, and "H, T-L and £ may depend on p^. We 
dehne 


A^' A'^ 


p'^ = (p°,o,o,p^), e = 


p- 


llz 


p 


Iz 


A" 


p” 


+ 2p^ ’ 


( 4 ) 


and t is the same as in the light cone GPDs, since it is Lorentz invariant. ^ dehned here 
approaches ^ in the light cone GPDs when hadron’s longitudinal momentum approaches 
inhnity. 

In the following we will focus on the generalized quark distributions in the nonsinglet 
case and consider quarks as external states. The on-shell conditions for the initial and hnal 
state quark 



= m 


2 


3 






and the definition t = lead to 


A° 

A^ 


= \ + pj — 


t 

4’ 


2ep 


2 

2 


m? + p1 


\J—At (1 — pl + t"^ — {^i'^pl + t) 

^4 {m? +p1 ) —t 


( 5 ) 


where we have kept a qnark mass m to regularize potential collinear singularities and have 
chosen A^ to point in the positive x direction without loss of generality. For |Ax| to be 
real, we have the following constraint 


—4t {l — + + f) > 0 



( 6 ) 


In the infinite momentum limit —?■ oo, this reduces to 

^ \/ —f + Am? ’ 

which is the constraint for ^ in the light cone GPDs. We will also assume ^ > 0. 


III. ONE-LOOP RESULT FOR GPDS 

In this section, we present the one-loop results for the quasi and light cone GPDs. As 
in the case of parton distributions, we choose the axial gauge = 0 throughout the 
computation since, in this gauge, the gauge link becomes unity. We also use a transverse 
momentum cutoff for regularizing the UV divergences. 

Let us start with the unpolarized case. From the definition of Eq. ([3]), it is easy to see 
that the quasi distributions yield the same result as the light cone ones at tree level 


{x, t) = {x, t, p") =6 (1 - x), 

{x, t) = {x, t, p^) =0. (8) 

At one-loop level, the contributing Feynman diagrams in the axial gauge are shown in 
Fig. [U Let us first look at the gluon-exchange diagram. In the axial gauge, the gluon 
propagator is given by —iD^y{k)/k'^ with the numerator 

(i\ 'fT'fj.ki/A -2 kfjky , , 

Df,^{k) = ^ - - -- + ra ^ 9 

n ■ k [n ■ ky 

where n ■ k = k^ = —1. The first term on the rhs. of the above equation leads to the 
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FIG. 1: One-loop diagrams for GPDs in the axial gauge. 


Feynman gange resnlt for the diagram, which can be written as 
d^k 




(27r) 


= -ig'^Ci, 


■u{p'‘ 

d^k 


-ig-f 


+ f — m 


7 


4 


-igY 


~'^g^iv 


— m 


ij) - ky 


tU 


u{p'‘ 


2Y 


[{k -h f - m2](p - /c)2 


8mk^ — — Ak^^ — 

[{k + 4)^ — m2][(/c — Ay — 777,2j(^p _ }^y 


u 






( 10 ) 


After a Feynman parametrization and integration over kP and /cx, we have the following 
resnlt for the hrst term in the cnrly brackets above 


hii 




y)p^ + (1 - y)^Y + (1 - vY'rrP 


( 11 ) 


where we have used k^ = xp^, and y is the Feynman parameter. 

The contribution of the second term in the curly brackets in Eq. flTOl) can be computed 
analogously, and the result is 


ri2 = 


g^Cp 

lOvr^ 


dx p^ / dy 


" 1 - 2 / 


dz 


0 [{k^ — {1 — z — y)p^ + {z — y)^y + Ayzp"^ + {z — yym'^] 2 


X — 2(1 — z — y)m{p^ + k^ — my^) — (1 — y)'y^A‘^ 

- {2k‘ - A')|(l - 2 - sjpV -(z- y)^7' - k‘l‘]}s(x - A. 


( 12 ) 


Performing the Feynman parameter integrations in Eqs. (ITTll and (IT^ . we are able to extract 
the contribution of Fi to the quasi GPDs "H and £ with the help of the Gordon identity. 
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The result reads 


Hi = 



ctsCp 

2 tt 


21n(m)(^+x) I ln(p^)($+x) ln(-t)(—2g+a:—l)(5+x) , (x—l)ln(^+l) 

(?+l)(x-l) 5K+1) 25($+1 )(x-1) 52-1 

ln(2)(c(x2+l)-2x(^+x-l)) ^ ?xln(l-|) ^ (l-x)xln(^) , 2Cln(|^) 
S($2-l)(x-l) ^ (C^-l)(x-l) + «K2_l)(x-l) ?^-l 

€ln(€+x) , «(«^-^)ln(f^|)-«2ln(^) 

«"-l ^ S($2-l)(x-l) 


41n(m)(x-^2^ 2(x-l)ln(p^) ln(-i)(- 2 C 2 +a; 2 +i) 

(C-l)(«+l)(^-l) «=-l (?2-l)(a;-l) 

I ln(g+l)(g+a:-2)(g+x) , 1 f _2_, (x-g)(g+x) \ , 

+ (g2-l)(x-l) + 2 l^g+1 (g2_l)(x-l) J 


+ 


(-2g+x-l)ln(-i^) ^ (g^+x)ln(^-l) 


2(g2-l) 


2g(g2-l) 



x—1 


ln(l-g)(-2g+x-l) 

f-1 
4(g+£)' 


X—1 


ln(_i«±2) 

2(x-l) 


X <-^ 


<X 


^ < X < 1 


ln(|^)-g(2g-x+l)ln(f4)+g(2g+x-l)ln(|^) 

2g(g2-l) 


X > 1, 


(13) 


Si {x,^,t,f^,Pl = 


asCp 'rri^ 

2vr {l-e)t 


0 

X < — ^ 


2K-l)(x-01n(,^)+2K2 + x)l,i|g 

+2^ (x + 1) In 4 ^ 2 (i _,^)2 

< X < e 

(14) 

4 (e" + x) In + 4^ (1 + x) In 

■^ < X < 1 


0 

X > 1, 



where /i is the transverse momentum cutoff for regularizing potential UV divergences. Some 
remarks on the above results are in order. To obtain Eqs. flT3|) and flTT)) . we take the limit 
/i —)■ cxo, and then —)■ oo, m —)■ 0, where we keep the leading term in p, and m of Hi and 
Si (for Si the leading term is (T(m^)), and ignore the power suppressed corrections of the 
type (1 /p^)'^{n >1). If we do not take p^ —)■ oo, power corrections should be kept in the light 
cone GPDs as well, in order for the quasi and light cone GPDs to have the same IR behavior. 
As in the PDF case [H, 13||, the quasi GPD result Hi does not vanish in the full x range. 
However, its collinear singularities exist only in the Dokshitzer-Gribov-Lipatov-Altarelli- 
Paris (DGLAP) and Efremov-Radyushkin-Brodsky-Lepage (ERBL) regions ^ < x < 1 and 
— ^ < X < ^. In the DGLAP region, the collinear singularities do not straightforwardly 
reduce to the corresponding PDF result with the Altarelli-Parisi kernel when taking the 
limit —?■ 0, because we assume a finite t in taking the limit m —?■ 0. If t = 0, as in the 

PDF case, the above assumption does not apply, and the term leading to ln(—t) above will 
lead to Inm^, and thus to the correct collinear behavior of PDF. S is UV convergent as 
expected, since it is zero at tree level. There is no UV divergence in H, but a logarithmic 
dependence on p^ instead, as in the one-loop results for PDFs. Moreover, the coefficient of 
Inp^ in H in the DGLAP region reduces to the corresponding PDF result when = 0. 
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The second numerator structure in Eq. ([9]) leads to 


To = C, 


d^k 


u{p”) (-*^ 7 " 




r 


4 ^ 




m 


X 


{p-k) n^ + n^ (p-k)^ 


(j)Z _ i^z'j _ j^y 


■u 


= -ig'^C, 


d^k 

Wy 


u{p") 


7 




7 "+ 7 ' 


1^ + ^ - m 


^ ^ — m 


r 


(j)Z _ i^z'j (j) _ i^y 


tU 




pz 


(15) 


and the result is 


■Ha (x,^,t,/i,p^) = 

( -^ -\ - ^ _In — H-_In — 


ctsCp 

271 




2(l+«) 


1 


+ In 


X — 1 

2(1-3;) ' (l-x)(l-$) x-^ 


2{x-e) 

{i-x){i-e) 


In^-i) + 


2 2 ^ (i-;i;)(i-«) 


(x-g) 4(x-g)(l-g)^ 


1 — X 


I i^+i) 4(x+g)(l+g)^ , 1 

■^(l-a:)(l+g) l-x n=p 

_1_1^ x-1 _ 3;+g ]„ x-i 

^ l—x (1—x)(l-g) X—g (1—x)(l+g) x+g 


X <-i 


-i<x<i 


^ <x <l 

X > 1, 


E2{x,i,t,p,p^) = 0, 


(16) 


where £ 2 {x,^,t, p,p^) is of order O^rrY/pl), which is power suppressed compared to £i in 
Eq. (ITT)) and therefore ignored. 

The last numerator structure yields 


T3 = Cf 

- 2 , 


= ig^Cp 


d^k 

WY 

r d^k 


u{p'') {-igY 


-u{p")- 


^ + ^ - m 


7 


. ^^^{p-k)^{p-k)^ k\ 

H. \ - 7zuip)d{x - 

\p-kf{p^-kY ’ ^ 


-f-m 


pz 
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j (27rY ' {p ~ k^Y {P^ ~ YY 

which contributes to "H {x,Yt,p,pY only with 

OsCf + pI(1 - P’' - |1 


(17) 




(18) 


27r (1 — xY 

Summing over all of these contributions, we obtain the following results for the gluon- 
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exchange diagram in Fig. [T] 




' ie+x)ln^^^ ^ {-2e+x^+l)ln^§^ ^ 


o^sCp 

271 


m^-i) 


2 ie-i)ix-i) 


P^(l — Xp 


In 5 + ihroSiT '■'|2(« + 1)] - iffhl ‘n 4^ 


a:+g 


(1_^)(1_^2) 




(1+C)(x-1) 

•2_9C2 


+ 


p(T-^ 


In "J + 2(ih)ri-V) (ln[16(^^ - 5^)1 - 2 In 

X+^^ 1 X^ _ 2(x-g^) I M 

2^(1—^^) 3^+^ (1—3:)(1—(1—3::)^ 

(«2+x)ln|^ (-2C2+x2+l)ln^ 


2«(€2-l) 2K2 _i)(^_i) 




+ 


p^(l—xP 


X <-^ 

-e < X < e 

^ < X < 1 

X > 1, 


(19) 


Now we present the one-loop results for light cone GPDs. As in the PDF case [l2j, the 
one-loop corrections for light cone GPDs can be obtained by first integrating over k^, then 
taking the limit —?■ cx) and integrating over k±. This leads to the following results for the 
three numerator structures in Eq. (jOj) 


Hi {x,^,t,p) = 


o^sCp 

27T 


m f Mi \ j_Eli_In f 

2?(l+«) \-tJ ^ (l-D(l+«) I -i 
x+e 1^ i+x^-2e 


1—x— 


^In 


1—X 


■ In C2l^2 + 


In — 


^25(1-52) inn ^2_3,2 2 {l-x){l-^'^) """ x-l-^ 

In ( H_J ( 

(1-52) m + (1-X)(1-C2) in 

l-x-2^ 1 1-x _ l-3:+2g 1 1-x 

1-^2 in in 


<x 

^ < X <1 
otherwise, 


r2(x-0i„/-n , 2C(i+x),„ Ae(i-x) 
-T+T (i+ 5 ) 2 ( 52 _^ 2 ) 


El (a;,^,f,p) = 


asCp Tn? 
271 —t 


i+« 

2 '' --L<i2 


1-f^ 

2 


1-^2 




^ < X < 1 
otherwise, 


H 2 (x, t, p) = 

X+^ 


oisCp 

2n 


(l-x)(l+0 

2{x-e) 


(1—a:)(l —m 


In -^ -|- In 

2 
2 


mr 




31+g 


x-^ 


(l-x)(l-0 


In 


{l-x){l+0 
1-0 2 


In 


(S 


l-xJ 


E 2 {x,^,t,p) = 0, 
H 3 (x,^,t,p) = 0, 


-^<x<^ 

^ < X < 1 
otherwise. 


( 20 ) 





























































the sum of which is 




asCp 

271 


n _L JL 

2?(l+€) 

_ l+a:^-2g^ Q 


n _L In 

-t + 25 (l-^ 2 ) 
(1-^)^ _ In 


In 


4g2 


52-x2 


(1+^)2 

l+x2-2g2 (l-x)2 


l+x2-2g2 1 /£^_ 

(l-x)(l-^2) -t (l_a;)(l-^2) 

0 


(1-^2) 


-^<x 
^ < X < 1 
otherwise, 




( 21 ) 


The results do not vanish only in the DGLAP and ERBL regions. There is a logarithmic 
UV divergence In /i^ in the above results, whose coefficient agrees with the coefficient of Inp^ 
in and also with the evolution kernel of nonsinglet GPDs in the DGLAP and ERBL 
regions. It is interesting to see that and are equal. This means no matching is 
required for the GPD E up to one-loop and leading power accuracy; therefore, the light cone 
E can be smoothly approached by the quasi S in the large momentum limit. The reason 
behind this is simple: the light cone GPD E and its quasi counterpart are zero at tree level 
and thus they are UV convergent at one-loop level. 

Next, we look at the contribution of self-energy diagrams in Fig. [1] The computation 
of the quark wave function renormalization factor is essentially the same as in the PDF 


case 


3 qu 

12| : the only difference is in the momenta of the quarks. Since the incoming and 
outgoing quarks now have different momenta, we have two wave function renormalization 
factors (p ± ^A) (or ^SZp (p ± ^A) for light cone GPDs). For simplicity, we denote 

6Zp(^p±^A) {6Zp {^p± &s 6Zp^±{^,t) {5Zp^±(^,t)). These factors can be related to 
each other as 


6Zp^p (e, t) = 6Zp,_ (-e, t) , 6Zp^p (e, t) = 6Zp,_ (-e, t) . ( 22 ) 

Using the same strategy as in Ref. [l2|, we obtain the following results for the wave function 
renormalization factors 




6 Zp2,- 


ctsCp 

277 


(y^sCp 

277 


x—1 

(1-S)- 

1—x 


In - A_ 

x-l 1-5 


dx < 




In 




+ 


In 


1 — X 
(1-0^ 
2G-5)2 _ 2(l-x) 

(i-x)(i-5)" (1-5)" 

In I 

Z77 in ^ 


4(i-5)T^-g) 

1—X 


(i-«) 

/ 2(x-5) 

(x-i)(i-5) 

2(^-5) 


In _l 3-2x-5 
x-l (a:-l)(5-: 


dx < 


(l-a;)(l-5) 


In 


(x-i)(5-i) 

I 1j^ 4(i-5)"(x-5) 


1 — X 


^-Ax+3 

(x-im-i) 

2(^-5) m _ 3-2x-5 

(x-i)(i-5) x-l (x-m-i) 


X 

^ < X <1 

X > 1, 

X < ^ 

■^ < X < 1 

X > 1, 


dZp^_ — 


cisCp 

277 


dx 


H a:+5—2 

(i-x)V 

n _ 1 _ 3^+5—2 

(i-x)V 

H _x+5-2_ 

(1-3;)V 


X < ^ 

^ < X <1 
X > 1. 


(23) 
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Summing over all these contributions and including also 6Zf,+i we have 


( 1 ) _ asCp 

^ ~ 2ti 


{m, I/) In fEf + /(-e, I/) In f^) - 

-/(-e, I/) In ^ - /(e, 2/) In iEf + /(-e, 2/) In 4(i+g)~4+^) 
+4/(-e,l/)-dF + T^ + ^ 
dy -(/(C, 2 /) + /(-^, y)) In ^ + /(^, y) In 4 (^_g^(^_g )2 

+/(-^)2/) In 4 (j^+^)(4+^)2 + 4:{f{^,y) + 2/)) + 


}ll y) In ^ - /(-e, 2 /) In g + 4 ^ + 


p 2 (l-y )2 


y <-i 

-i<y<i 


i<y<l 

2 / > 1 , 


where 


/(^, 2 /) = 


2 / 


1 -e 1 - 2 / 2 ( 1 - 02 ' 

Similarly, on the light cone one obtains 


(24) 


SZpi — 


6 Zf2 ,- — — 


(^sCf 

27r 

(^sCf 

27r 


1—a: 




2 (^-g) 
dx 


T (In - 2/-2 + 2 In 1 —^ 

^ \ ' 2 _ X 


2(^-0 

(l-x)(l-5) 


^ < X < 1 

otherwise, 


In ■ 


+ 21 nl 3 |j ^<x<l 

otherwise. 


SZfs- — 0 , 

and the complete result including SZf,+ is 


(25) 


7 ( 1 ) _ 


OsCf 

2f 


dy < 


-/(-{,9)ln^ - 2/(-J,!/)lni^ + ^5“!^ -( < V < i 

-m, y) + n-i ,»)) In ^ - 2 /( 5 . y) 111 
-2/(-5.9)l,i^ + .iij + ^-jlj 5<9<1 

0 otherwise. 


From the above results, it is clear that the coefficients of In ^ in ^ and In ^ in zl!^ 
in the DGLAP region [0 1] reduce to the PDF result when 0. Several further checks 
have been performed on our results: the one-loop results for GPDs reduce to that for 
PDF when 0^ ^ 0 (although for H, this is non-trivial); the x integrals of the Feynman 
part Eqs. (fTOj) contribution to GPDs, which is the most complicated one among all three 
numerator structures, agree with the form factors of a local quark current computed in the 
Feynman gauge. We also checked the polynomiality of the logarithmic terms in the above 
results. 

In the case of polarized GPDs, the vertex contribution can be obtained by replacing 7 ^ 
with 7 ^ 7 ® in Eq. flTU]) . whereas the self-energy contribution remains the same. We present 
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the result of the vertex contribution below. For the quasi GPDs, we have 




agCp rn? 
2 vr 


0 

X < —^ 

2(5-1) (5-i^)ln(;^)+2 e(l+i)lng 

,0 55-2 , M (l+«)^(«"-^^) 

+2 (5"+ 3,) In 1 

-i<x<i 

4e(l + x)ln (^^)+4(e2 + x)lni^ 

i <x <l 

0 

X > 1, 


(26) 


whereas on the light cone we have 

H^^\x,^,t,p,p^) = H^^\x,^,t,p,p^), 

(27) 


IV. ONE-LOOP FACTORIZATION 


Before we construct the factorization formula connecting the quasi and light cone GPDs, 
let us summarize the one-loop results in the previous section. For the unpolarized case, we 
have 

n{x,^,t,p,p^) = [1 + ]^{5Zf{p - ]^/\) + 5Zf{p + ]^/^))]5{l -x)+H^^\x,^,t,p,p^) 

Ki) 


= {I + Z^p’)5{l - x) 

H{x,^,t,p) = [1 + ^{5Zf{p- ^A) + 5Zf{p + ^£^))]5{1 - x) + H^^\x, p) 
= (1 -h Zp'^)6{l - x) + H^^\x,^,t, p). 


(28) 


Similar results can be written down for the polarized case. A crucial difference between "H 
and H is that the latter vanishes in the regions x < and x > 1, whereas the former does 
not. The connection between the quasi and light cone GPDs can be established as 




(29) 


up to power corrections suppressed by p^, where the integration range is given by the support 
property of the light cone GPD. 

The matching factor Zh can be perturbatively expanded as 


.y y p 


«5 y{l) (x i P 


'-Z)^> ( 21 1 + h.o., (30) 

where h.o. denotes higher-order contributions. From the one-loop results for PL and H, the 
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matching factor can be extracted as 


Zh\vX,p^Ip‘)ICf 


= 


(C2+r,)ln^ , (-2C2+r?=+l)lni^ ^ ^ 

2C(C2-1) + 2(C2-l)(r,-l) + p-(l-r,)2 

1_l^2_9^2 


(1 + h + ( M4(1 - .()=] - In 

^+c 


2C(1+C) 

^+c- 


^+c \A(i^2 _ 2\] J _2-K_I_ I 

2c(i-cb ^ h JJ + (i+c)(»7-i) ^ Phi 


l+rj2-2C4 


-r ;)2 


(i_^)(i_C^ In + 2 (l-;)(i-V) (- (2)] + 2 ln(l - r;)) 

V+C^ In It< — 2(??-C^) I ti 

2C(i-Cb p+C (i-piCi-C^) pH^-v)^ 


(C"+P)ln^ (-2C"+r?2+l)lni^ 


2C(C^-1) 


2(C^-1)(P-1) ^ P^(l-»7)^ 


h < -c 


-C <v<C 


C < h < 1 

?7 > 1. 

(31) 


The above matching factor is valid for y > ^, however, it can be extended to the full y range 
as 




M 


y y p J 2/L \y y p J 

+ F2 < ^ < < y) 

\y y p"J 


X i y 


y y p 


+ ^3 < a: < I/) + F4 — e{x > CjO{x > y) 


X i y 


y y p 


(32) 


where -Fi, 2 , 3,4 are given by the matching factor in the four different regions in Eq. (13T|) ( with 

a replacement Ina —?■ l/ 21 na^ so that they are real functions); shall be replaced by yP^ 

with being the averaged longitudinal momentum of the external hadrons. The validity 

2 

of the above equation can be checked by explicit computations. The coefficient of In p- is 
the same as the evolution kernel of the light cone GPD. 

Near p = I, one has an extra contribution from the self-energy correction 

4 " {-,-,-}]=6Z^ff(2xlas)i(l- v) (33) 

\y y p"J 

with 


6Z 


( 1 ) _ 
H ~ 


c^sCp 

2ti 


dp < 


(/(C,h) ln^ + /(-C,^) In^f) - 

-/(-C, h) In jl - /(C, V) In ^ - /(-C, P) ln[4(C + ??)(!- p)] 

+4/(-C,h) - - t4 + t4 + p4i%)2 

-(/(C, v) + /(-C, v)) In4 - /(Cl v) - 0(1 - 0] 

-/(-C, 0 ln[4(C + 0(1 - 0] + 4(/(C, 0 + /(-C, 0) + 

-}(C, 0 iV^ - /(-C, 0 In + ^ + ^ 


V <-C 


-C <v<C 


C <p <1 
p > 1 . 

(34) 


One can check to see that <5^4 provides a plus distribution for the singularity at x = y 
in the matching factor Z^j^\x/y,^/y, p/p^). The above matching factor also transforms the 
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logarithmic dependence on in t,/i,p^) into the renormalization scale dependence 

in H (x, t, p). Note that we have the same linear divergence as in the PDF case. Moreover, 
when ^ 0, the region disappears and the matching factors in the remaining regions 

reduce to those for the PDF. 

In the above result, we take into account the quark contribution only. The antiquark 
contribution is given by making the following replacement in Eq. fl32p 

X —X, y —y. (35) 

Summing over both the quark and the antiquark contribution, one obtains the complete 

2 

matching factor, which contains in the In % term the complete evolution kernel of the light 
cone GPD. 

The factorization of E can be constructed analogously. However, from the results of £ 
and E, the matching factor for E is simply given by 

ZE{x/y) = 6{x/y - 1) (36) 

up to one-loop order and leading p^ accuracy. Since E does not show up at tree level, it 
is UV convergent and thus does not have a cutoff dependence. Accordingly, £ does not 
have a logarithmic dependence on p^. Therefore, the light cone GPD E can be smoothly 
approached by the large momentum limit of E. This is also true, in principle, for other light 
cone quantities that do not exhibit a UV divergence. The simulation of such quantities on 
the lattice are, therefore, relatively simple. 

Since = E[^^\ and E^^^ = £^^\ the above factorization also applies to 

the polarized GPDs with the same matching factors. 

In the following, we consider the matching for the distribution amplitude, whose evolution 
kernel can be obtained from that of the GPD as a limiting case. Here we focus on the simplest 
type, the distribution amplitude of the pion. The light cone pion distribution amplitude 0(x) 
is given by 


where the two quark helds are separated along the light cone, and x (1 — x) denotes the 
momentum fraction of the quark (antiquark). As in the case of GPDs, it can be studied 
from the large momentum limit of the following quasi correlation 


with the two quark helds separated along the spatial direction. The one-loop factorization 
for the pion distribution amplitude can be written down analogously as 





dy Z^{x,y,p,p^)(j){y,p). 


(39) 


The matching factor Z^{x,y, p,p^) can be obtained by starting with the light cone and 
quasi distribution amplitudes, Eqs. flTT)) and fl38|) . and computing the one-loop corrections, 
respectively. It can also be obtained from the above matching factor for GPDs by crossing the 
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initial quark to the final state, which leads to the following replacement for the momentum 
fractions 

C ^ v/C ^ 2 X- 1 . (40) 

2y-l 

This corresponds to setting 

x^2x-l, y ^2y-l, ^ ^ 1 (41) 


in Eq. (132|) . Note that one also needs to replace the averaged longitudinal momentum of 
external hadrons below Eq. (13^ by p^/2 with being the longitudinal momentum of 
the pion. 

Expanding the matching factor Z^{x,y,pi,p^) as 


Z4x,y,ii,p’) = S(x -y)+ y,j/) + .... 


(42) 


e then have 
hi) 


z), ’{x,y,fi,p'^)/CF = Gi(x,|/,p,p^) 6 '(a: < 0) + (^ 2 ( 0 :, |/, p,p^)6'(0 < x < y) 

+ Gsix, y, /i, p^) 6 »(?/ < X < 1) + G 4 {x,y,fi, p^)9{x > 1) (43) 


with 

^ , X 1 —Xx, X , X 1 —X 1 X , (x — 1 /)^ 

Gi x,|/,/i,p = (—-- + ——) In-- + (—-- - —— +-) In—-— 

^2{l—y) 2y ' x — 1 ^2{l — y) 2y y — x x(x — 1) 


p^{x — yY ’ 


^ /X — 1 1 X - / X 1 —X 1 r./ 421 

G,(x, 9 ,^,P ) = (_ + _) In - + 1„|4(X - y) ] 

i_ + £^) ,„[4:.(1 - :.)] + ( " ' + _ ^) l„i^ 

' ^2{y — l) 2y y — x 


2{y - 1 ) 2y 

11 a 

+ - + ^ 


y y — x p^{x — yY' 

G3{x,y,p,pY = ^ 2(1 - X, 1 - y,y,pY, 

, X 1-XxX-l , X 1-x lx, x(x -1) 

G,(x,y,ix,p ) = __ 


+ 


p^(x — yY ’ 


(44) 


Here we have taken into account charge conjugation invariance and the fact that one-loop 
diagram for the pion simultaneously involves the quarks and the antiquarks state. We have 
checked to see the above matching factor agrees with the result of direct computation. 

Near x = y, one has an extra contribution from the wave function renormalization 


Z^^\x,y,y,pY/GF = 5Z^^\2 ti/ as)5{x - y), 


(45) 
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where is given by 


= 


c^sCf 

27r 


+ 


' 7 ’_ 1 

'2{1-J/)" 

1 __ 

2(y-l) 2y p^(i/-a;)2 


2 (?/-l) 


In + (; 5 ^ + f - ^) In ^ 

V—X' y—X ' iu V — X' x — v 


y-x/ 

y- 


' 2 p' 


2 y 


y-xJ 




1 


dx < 


+ ^)ln^ 

y—X'' x—\ 


2y ■ y-x 

-y 1 


2(l-p)2 2(p-l) 


1—x 


_L ^ I J_ 

2(j/-l) ^ p2 -T 2y 




+ 


2(t=p)2 + + + 

I 1 I I 1 I 1 I M 

2?/ (1-y)^ 2(l-y) y-x ^ p^{y-x)‘^ 


+ -In 

2y y-x/ 


x—1 


x-y 

X 

1 


1 


+ —t —tT-) “ 

j/-_i \ y—X' y—x 

a 


2y‘ 




y-x/ 


In ■ 


x-y 


+ 1 .+ 


2(!/-1) 2p p^{y-x)‘/ 


X < 0 
0 < X < y 
y < X < 1 

X > 1. 

(46) 


One can explicitly check to see that the wave fnnction renormalization factor provides a pins 
prescription for the factor in Eq. fl43l) . 


V. CONCLUSION 

We have presented the one-loop matching conditions for the nnpolarized and polarized 
generalized qnark distribntion in the nonsinglet case. The matching conditions relate the 
qnasi GPDs dehned in terms of spacelike correlations and the light cone GPDs. For the GPD 
H (H), the matching is constrncted in analogy with the PDF matching, and the matching 
factor rednces to that of the PDF in the limit 0. For E (E), as it is UV convergent, the 
matching factor is trivially given by a 5 fnnction, implying that the light cone GPD E [E) 
can be smoothly approached by its qnasi connterpart S in the large momentnm limit. 
This facilitates its extraction from lattice simnlations. We have also presented the matching 
condition for the pion distribntion amplitnde. 
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